Quantum Theory of Cavity- Assisted Sideband Cooling of Mechanical Motion 
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We present a fully quantum theory describing the cooling of a cantilever coupled via radiation 
pressure to an illuminated optical cavity. Applying the quantum noise approach to the fluctuations of 
the radiation pressure force, we derive the opto-mechanical cooling rate and the minimum achievable 
phonon number. We find that reaching the quantum limit of arbitrarily small phonon numbers 
requires going into the good cavity (resolved phonon sideband) regime where the cavity linewidth 
is much smaller than the mechanical frequency and the corresponding cavity detuning. This is in 
contrast to the common assumption that the mechanical frequency and the cavity detuning should 
be comparable to the cavity damping. 

PACS numbers: 42.50.Lc, 07.10.Cm, 42.79.Gn 



Recent years have seen significant advances in fabri- 
cating, measuring, and controlling mechanical systems 
on the micro- and nanometer scale [H, 0, 0, 0, 0, Sj- 
These advances have opened the possibility of observ- 
ing quantum effects in mechanical devices. A num- 
ber of experiments to explore the quantum regime of 
these systems have been proposed, including measur- 
ing their coherence and entangling them with other sys- 
0, 0, US U, 12 1, generating nonclassical states 



terns 

13 , Il4j , and making quantum- limited measurements 
15l . [16[ . Many of these proposals require the mechanical 



system to be cooled to its ground state, i.e. to tempera- 
tures below 20 mK even for 1 GHz resonators. This is dif- 
ficult or impossible using bulk refrigeration, and a more 
promising means is to make use of non-equilibrium cool- 
ing techniques, analogous to the laser-coolin g sc hemes 
developed for trapped ions and neutral atoms [171 ] . 

Such schemes fall into two categories. In the first, cool- 
ing is achieved via an active feedback loop which is used 
to cancel the cantilever's thermal motion 18. Hp , I20I ]. 
This approach (including its quantum limits [211 ] ) has 
been considered elsewhere. Here we will focus on the 
second category: passive, non-feedback-based cooling. 
In this approach the cantilever displacement is coupled 
paramctrically to a driven resonator (or two-level sys- 
tem). When the frequency of the drive applied to the 
resonator is chosen appropriately, this parametric cou- 
pling ensures that the cantilever is preferentially driven 
to lower energy states. The lowest temperature to which 
this process can cool the cantilever is determined by the 
resonator's quantum fluctuations (photon shot-noise), 
which drive the cantilever randomly and compete with 
the cooling. 

This balance has been considered theoretically for a 
few specific realizations of the cooling system: a Cooper- 
pair b ox 12 211 . the superconducting single-electron transis- 
tor [S 123LI24I] , quantum dots [HI and ions [26| . In this pa- 



per we present a quantum mechanical description of the 
self-cooling of a cantilever coupled to an optical cavity via 
radiation pressure. This system is of particular interest 
because of its relative simplicity - the cantilever is cooled 
by a single mode of the electromagnetic field. It is also 
of interest because this system has already been realized 
experimentally, and recently produced very promising re- 
sults on self-cooling using both photothermal forces and 
radiation pressure forces 

di US US US EE US ES E3 ■ 

Cooling rates and steady state temperatures for such 
a scheme (and of related systems like a driven LC cir- 
cuit [351 ]) have been derived only (semi-)classically so far. 
Most notably, Braginsky and Vyatchanin [3(| considered 
the optomechanical damping within a scmiclassical the- 
ory. Here, we use the quantum noise approach to find 
simple and transparent, but fully quantum-mechanical 
expressions, starting directly from the spectrum of the 
force fluctuations. Our results are valid both for the 
good-cavity regime (resolved mechanical sidebands) and 
the bad-cavity regime (unresolved sidebands). Earlier 
(semi-) classical calculations [2^, EH, EH are reproduced 
in the appropriate limits. We show in particular that 
it should be possible to cool the cantilever to its quan- 
tum mechanical ground state by choosing the cantilever 
resonance frequency much larger than the cavity ring- 
down rate, a regime that has not been considered so 
far. Cooling the cantilever to its ground state will en- 
able the realization of many of the aforementioned tasks, 
as well as investigations into the optomechanical insta- 
bility 



37 



38|, [39( in the quantum regime. The results arc 
complemented by an exact solution of the coupled equa- 
tions of motion to account for the "strong cooling" limit, 
where the cooling rate exceeds the cavity ringdown rate. 

Consider a mechanical degree of freedom x coupled 
paramctrically with strength A to the cavity oscillator 



H = h(uj R -Ax) [a^a~{a 1 a)]+H M +H dlivc +H K +H r (1) 



'At 



2 



where is the cavity resonance frequency at the equilib- 
rium cantilever position x — in the presence of the mean 
radiation pressure, Hm = huiM^c is the mechanical os- 
cillator, -f/drive is the optical drive, H K represents the cav- 
ity damping, and Hr represents the mechanical damping. 
For a cavity of length L, the change in cavity length with 
x gives a radiation pressure coupling A = H-cvr/L. 

Opening the port used to supply the classical drive to 
the cavity also admits vacuum noise. Writing the cavity 
field in terms of its classical and quantum parts a(t) = 
e _ta 'i'*[o+d(t)] yields the photon number autocorrelation 
function 



Snn(t) = (fit(t)6(t)fit(0)6(0)> - <fit(t)6(t)) s 

a iAt-f\t\ 



ne 



(2) 



where A = lj-l — lor is the laser-cavity detuning and n is 
the mean photon number. Because the number fluctua- 
tions are the result of interference between the classical 
drive amplitude and the vacuum fluctuation amplitude, 
they decay at the amplitude decay rate k/2, not the en- 
ergy decay rate k as is sometimes naively assumed. 
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FIG. 1: (Color online) (a) Noise spectrum of the photon num- 
ber in a driven cavity as a function of frequency when the 
cavity drive frequency is detuned from the cavity resonance 
by A = +5k (leading to heating: dashed line) and A = — 5k 
(cooling: solid line), (b) Effective noise temperature T e s as 
a function of the detuning of the drive from the cavity reso- 
nance, see Eq. ((5j), for a cantilever frequency matching the 
detuning: lom = |A|. The noise temperature is positive on 
the left side of the resonance, where optomechanical cooling 
produces a minimum reachable cantilever phonon number 



from the cantilever, while noise at negative frequency 
5ff(— wm) corresponds to the ability of the cavity to 
emit a quantum of energy into the cantilever. Assuming 
lum 3> r^/, the net optical damping rate of the cantilever 
from Fermi's Golden Rule is 



Topi — T2 [Sff(^m) — Sff(~ Wm)] £zp] 



(4) 



where xzpf is the cantilever zero point position uncer- 
tainty, and Sff = fi 2 A 2 S nn . 

As can be seen from Eq. ^ and in Fig. (QJi) , for posi- 
tive detuning, the noise peaks at negative lo meaning that 
the noise tends to heat the degree of freedom x. For nega- 
tive detuning the noise peaks at positive lo corresponding 
to the cavity absorbing energy from x. Essentially, the 
interaction with x (three wave mixing) tries to Raman 
scatter the drive photons into the high density of states 
at the cavity frequency. If this is uphill in energy, then 
x is cooled. A somewhat similar mechanism involving 
an oscillator with non-linear damping was considered by 
Dykman 0. 

The radiation pressure noise is completely non- 
equilibrium but can be assigned a unique effective tem- 
perature assuming that the mechanical oscillator is per- 
fectly harmonic with sufficiently weak optical and me- 
chanical damping, Tm + T opt <C k, so that it responds 
to the noise only at frequencies ±Wm- In the absence 
of mechanical damping, the mean number of mechanical 
quanta present in the steady state of the cantilever 
is given by the detailed balance expression (for A < 0) 

+ 1 _ S F f{+um) 



M 



Sff(— wm) 



exp 



M 



(5) 



The noise temperature T e g will become negative when 
the noise leads to heating (see Fig. []}. From Eq. ([3]) we 
obtain 



' M 



(com + A) 2 + ( K /2f 
4jmA 



(6) 



(black line) corresponding to T e g, see Eqs. (|5I7 



For the special case of detuning A = — lom, we have the 
simple limit 



The power spectrum of the noise is 



4w M 



(7) 



+oc 



dt, 



Snn{t) — n 



(w + A) 2 + (k/2) 2 



(3) 

The symmetric part of the noise spectrum was discussed 
in [361 ] to estimate the additional noise introduced by 
the optomechanical damping mechanism. Note however 
that it is the asymmetry in the noise which leads to cool- 
ing or heating. The noise of the radiation pressure force 
F = hAh at positive frequency, Sff{+^m), corresponds 
to the ability of the cavity to absorb a quantum of energy 



which shows that in the limit com S> «, the quantum 
ground state can be approached provided that the optical 
damping dominates over the mechanical damping. For 
the case A = —u>m, we obtain from Eq. (|4]): 



r op t = 4 



V L 



1 



1 + f-S-V 

' V 4u M / 



(8) 



If the mechanical damping Tm is not negligible compared 
to r opt , then a rate equation yields the full expression for 
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the mean number of mechanical quanta in steady state 
(Fig. ®, 



n M 



r 



opt^-M 



r M n M 



■ opt 



M 



(9) 



where is the equilibrium mechanical mode occupation 
number determined by the mechanical bath temperature. 
In the limit of large mode occupation we can replace the 
occupation numbers by the corresponding temperatures 
and we see that this expression reduces to the expected 
classical expression for the final effective temperature. 

We emphasize once more that using a large value of 
the detuning |A| = wm 3> k offers the advantage of be- 
ing able to reach an arbitrarily small minimum phonon 
number by optomechanical cooling, sec Eq. ([7]). The only 
price to pay is to increase the input intensity, in order to 
keep constant the number of photons inside the cavity, 
n = n max /(l + (2A/k) 2 ), which determines the cooling 
rate. Here n max is the photon number at resonance, pro- 
portional to the input power. (Note however that n is 
limited by our weak coupling assumption r op t+rAf <C k) 
In the good cavity limit, the optimum cooling condition 
A = u>m requires detuning the drive laser by more than 
a cavity linewidth. Although this decreases the cooling 
efficiency in terms of the input power, it increases the 
efficiency in terms of the circulating power. This is ad- 
vantageous in experiments where the limiting factor is 
heating of the cantilever by residual absorption from the 
circulating power. 

Apart from the cooling or heating of the cantilever, 
there is the well-known optomechanical frequency shift 
of the mechanical mode ("optical spring effect") (see e.g. 
34, H^), which can also be expressed in terms of the 
force spectrum, using second-order perturbation theory: 



5u> 



M 



J ZPF 



1 



1 



LUM 



(10) 

We emphasize that Eqs. ©, ©, ©, and (TO]) can 
be used to obtain the optical spring effect, the cooling 
rate, and the cooling limited phonon number for arbi- 
trary force noise spectra. Such more complicated spectra 
might, for example, result from the contribution of more 
than one mode, or relate to a model different from the 
one discussed here, e.g. the coupling of the cantilever to 
some electrical circuit [35|. A related discussion in the 
context of coolin g by a single-electron transistor has been 
provided in Ref. |23j . 

We now compare with commonly employed simpler 
models [27, 3(1 35|, where one postulates the cavity light 
intensity to approach its position-dependent equilibrium 
value in an exponential fashion, dn/dt = 7(77,(2:) — n) 1 
with n(x) denoting the equilibrium photon number as a 
function of position. Linearizing this equation and solv- 
ing for the motion of x yields an effective optomechanical 




2 3 
Detuning co m /k 

FIG. 2: (Color online) The steady state phonon number um 
obtained through cavity sideband cooling as a function of de- 
tuning, with A = —ujm. Thick lines represent the weak- 
coupling quantum noise result, Eq. (5), for three different 
laser powers (expressed via constant n max , so the actual pho- 
ton number n increases towards smaller u>m)- The bath tem- 
perature corresponds to n M — 100 phonons, and the op- 
tomechanical coupling has been fixed to yield a realistic ra- 
tio r op t(n = 1,um — » 00) /Ym — 0.1. The dotted red line 
shows n M = (k/(4o;a/)) 2 from Eq. J7J). The dashed (dash- 
dotted) lines show the results of the input-output equations, 
for Ym/k = 10" 4 (10" 5 ) and n max = 10 6 . They were obtained 
by integrating the exact spectrum S cc (uj), Eq. (|14[) , which is 
shown in the inset as a density plot, for Ym/k = 10~ 4 (Arrow 
indicates splitting of resonance; "unstable" regime does not 
permit A = —ujm)- The weak-coupling limit is recovered for 
Ym/k -> 0. 



damping rate 27, SUE 



r' 



hA 2 



On 



opt 



7717 1 + (wm/t) <9A 



(11) 



In general, this is not equal to the correct quantum- 
mechanical result @, if we use a fixed 7 oc k, and it 
would predict an optimum damping rate for 7 compara- 
ble to u>m- For damping by radiation pressure, the full 
analysis shows that we must instead employ a detuning- 
dependent effective decay rate, 7 = ((f) 2 + A 2 )/k. Then 
we recover Eq. ^ in the limit % -C n. We note 
that Eq. (TTT)) is always valid when applied to bolometric 
forces, where one describes the time-lag of the bolomet- 
ric force due to a finite heat relaxation rate 7 <C k, and 
where a quantum theory would be different from the one 
presented here due to the dissipative nature of the force. 

The quantum result jl]) for r op t can be reproduced by 
a classical theory based on linearizing the equation of 
motion for the complex light amplitude a itself (33l . |3(| , 
of the form da/dt = i(A + Ax)a — -|(a — a). Of course 
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it is still not possible to recover the correct formula ([9]) 
for the steady-state phonon number within such a purely 
classical theory, unless it is extended to include the zero- 
point fluctuations. 

In order to have direct access to the mechanical and 
optical fluctuation spectra, we now derive an exact solu- 
tion of the linearized Heisenberg equations of motion for 
d and c, where we use the input-output formalism [4l| to 
take into account the damping and the fluctuations, both 
for the optical and the mechanical degree of freedom (see 
also jHU) : 



d = iAd — — d — V^^in + ict(c + c^) 



-ILJmC 



M 



(12) 

i(a*d + <*«?). (13) 



Here the effective light amplitude has been expressed in 
terms of a frequency, a = a(u> rxzpf / L) with |a| 2 = n. 
The solution yields the cantilever spectrum S cc {lo) = 
J dte™* (ct(t)c): 



I I 2 

O / \ r M<Jth(w) + ^-cr pt(^) 
Scc{u) = rrr, ^ 2 , (14) 



WW 



where 



r o pt/re > 1/2, the mirror resonance peak splits into a 
pair of peaks at -wm ± a, getting hybridized with the 
driven cavity mode (see arrow in inset of Fig. [5]). At even 
larger photon number or smaller u>Ui when a ~ Wm/2, 
the static bistability [42[ precludes reaching the desired 
detuning A = -wj/ (see inset of Fig. [5]). Thus, we 
note that the far-detuned regime u>m 3> t has the addi- 
tional strong advantage of avoiding the bistability, which 
already interferes with cooling in some current schemes 

M- 

We have obtained the full quantum theory of cavity 
assisted sideband cooling of a cantilever, based on the 
quantum noise approach applied to the fluctuations of the 
radiation pressure. In the previously unexplored regime 
of a detuning much larger than the cavity linewidth, 
the cantilever may be cooled to arbitrarily small phonon 
numbers. The theory analyzed here may form the basis 
for the interpretation of future optomechanical experi- 
ments that will venture into the quantum regime of me- 
chanical motion. 
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othM = (n^ + l)|SM| 2 +^| x -i( w ) + iS(a;Xl5) 
cr opt H = k 2 \ X r(u)\ 2 \Xm(u)\ 2 (16) 
Af(w) = ^Mxi/'W+Mw). (17) 

We introduced the response functions of mirror and 
optical resonator, xm(&) = [^f 1 — *( w — ^^/)] 1 and 
Xr(u) = [f — i(u + A)] -1 , and we defined the optome- 
chanical "self-energy" £(u>) = — i\a\ 2 (xn(to) — Xr(~ u ))- 
The quantum noise results given previously are valid 
in the weak-coupling limit rj\/,r opt -C K- Then, 
the optomechanical damping and the "optical spring" 
frequency shift can be read off the self-energy as 
ImS(wj/) = — r opt /2 and RcS(o;jv/) = 5ojm, coincid- 
ing with the expressions given above (Eqs. ^ and 
(|10p ). The steady state average phonon number njyi = 
/ flF'-' cc ( a; ) reproduces Eq. ([9]). The optical output spec- 
trum (of d ut = d[ n + y/Rd) displays a Stokes peak at 
u> = —ujm and an anti-Stokes peak at w = +wm, with 
weights given by r op t(nM + l)n^ and T opt n M (nM + 1), 
respectively, which are the rates of processes leading to 
heating/cooling of the cantilever by red-/blue-shifting a 
reflected photon. 

The exact solution allows us also to discuss the regime 
of strong cooling. In the limit r opt « 4|a| 2 /K 3> Tm, it 
adds a term n^TM/ft + 2n^r opt /K to Eq. (J9j) for um 
(assuming k, a <C (jJm and Tm *C k). Then a min- 
imal phonon number n 
function of r opt at (r opt /re) 



min \ 
M — 
2 _ 



n^rjvi/re is found as a 
(r M /re)(nT/n°)/2. For 
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